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ABSTRACT 


Several  scaling  generalizations  for  a  pulsed  CO  FDL  are  presented 

and  discussed.  It  is  shown  how  all  necessary  parameters  and 

quantities  of  interest  can  be  scaled  to  analyze  a  specific  problem 

in  terms  of  more  general  classes  of  results.  The  sensitivity  of 

the  system  response  and  the  characterization  of  the  gain  saturation 

as  a  function  of  total  radiation  intensity  are  investigated.  It  is 

shown  that  the  transient  quantum  energy  and  power  conversion 

efficiencies  are,  to  a  good  approximation,  mainly  dependent  only 

upon  the  total  energy  deposition  as  a  function  of  time.  For  typical 

cases  of  interest,  the  scaled  cavity  loss  coefficient  is  small  enough 

that  oscillation  is  high  above  threshold,  and  it  is  found  that  turn -on 

occurs  (typically)  after  ~ 0.6-1.  0  J/liter/torr  has  been  deposited,  and 

that  steady  state  is  attained  after  —  1.  9-2.6  J/liter/torr  has  been 

*  \  j 

deposited.  (These  energies  increase  approximately  as  W  ,  where  We 

2  6 

We  /p cq)*  At  steady  state»  the  scaled  total  saturation  intensity 

is  ~  0.6-6.  4  W/cm  /torr  over  the  range  of  temperatures  from 
60-300°K. 
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1.0 


INTRODUCTION 


A  molecular  kinetic  model  for  an  electrically  excited  (CO,  N2,  Ac, 

He,  ,  ,  .  )  gas  mixture  has  been  constructed  for  the  calculation  of  the 
radiative  characteristics  of  a  CO  laser  oscillator  or  amplifier,  A 
computer  program  for  the  numerical  solution  of  the  steady  state  or 
transient  master  ecmadon  for  the  diatomic  species  has  been  developed 
to  predict  vibrational  population  distributions,  gain  and  saturation  para¬ 
meters,  energy  transfer  and  extraction  rates,  conversion  efficiencies, 
output  intensities,  and  spectral  distributions.  The  radiative  calculations 
are  self-consistent  with  the  saturated  gain  on  all  oscillating  transitions 
equal  to  the  losses.  The  plasma  characteristics,  described  by  a  Boltzma.m 
distribution,  are  adjusted  self-consistently  in  the  calculations  as  a  function 
of  time  to  account  for  electron  heating  for  a  given  temporal  input  power. 

A  more  complete  discussion  of  the  assumptions  upon  which  this  model  are 
based  has  been  given  elsewhere  ’  and  numerical  results  for  a  variety  of 
operating  conditions  were  presented  and  compared  with  experimental  data 
from  a  pulsed  high  pressure  E-beam  stabilized  electric  discharge  CO 
laser  oscillator  (EDL).1  Furthermore,  very  good  agreement  has  also 

been  obtained  in  comparison  of  theoretical  calculations  from  this  model 

3 

with  results  of  transient  gain  relaxation  measurements  and  with  small- 
signal  gain  measurements  over  a  large  number  of  vibrational  bands. 

It  is  possible  to  make  several  reasonable  approximations  that  allow 
certain  general  scaling  laws  for  a  CO  EDL  to  be  formulated.  Because 
of  the  successful  comparison  of  predictions  of  the  theoretical  model  with 
experimental  data  from  the  E-beam  pulsed  device,  the  results  of  scaling 
the  transient  analysis  for  an  electrically  excited  pulsed  CO  oscillator  will 
be  presented  here.  These  generalizations  depend  only  upon  a  few  basic 
scale  parameters,  which  make  it  possible  to  predict  the  temporal  output 
characterise -s  of  a  large  class  of  systems.  However,  the  present  results 
are  limited  to  gas  mixtures  containing  CO  and  monotomic  diluent  gases 


2 


only.  This  restriction  provides  a  considerable  simplification  insofar 
as  the  generalized  results  do  not  depend  upon  the  ratio  of  two  diatomic 
species  present;  moreover,  it  is  necessary  at  the  present  time  due  to 
uncertainties  in  the  (CO,  N2)  and  (N2,  R.)  cross-relaxation  rates  as 
a  function  of  kinetic  temperature  and  vibrational  quantum  numbers. 
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2.0 


MOLECULAR  KINETIC  MODEL:  SUMMARY 


The  physics  of  the  present  theoretical  model  that  is  used  to  predict 
the  characteristics  of  a  transient  pulsed  electrically  excited  CO 
laser  oscillator  is  described  by  four  basic  systems  of  equations  that 
couple  the  CO  vibrational  state  population  densities,  the  radiation 
intensities,  the  gas  kinetic  temperature  and  the  plasma  kinetic  para¬ 
meters  (electron  density  and  temperature)  as  a  function  of  time.  The 
general  formulation  of  these  equations  and  assumptions  made  will  be 
briefly  summarized  here,  and  it  will  be  shown  where  certain  approxima¬ 
tions  (most  of  which  are  not  required  for  the  complete  computer  analysis, 
however)  make  it  possible  to  obtain  essential  results  in  terms  of  a 
very  small  set  of  scale  parameters.  In  fact,  two  basic  parameters 
related  to  rumping  and  cavity  loss  are  all  that  will  be  required  to  present 
a  very  generalized  description  of  the  (CO,  X)  EDL  under  the  set  of  approxi¬ 
mations  to  be  discussed  below. 


The  physics  of  the  molecular  kinetic  processes  is  described  by  the  master 
equation  for  the  vibrational  population  densities  of  the  CO  molecule, 


Wdt  -  ne  I  [nv,  <  v(E)  vUv(E)  >  -  ”v  <v(E)  cr  v_v.  (E)  >  ]  + 
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,VT 


2  NX  v \  l-v(C°’  X)  "v  +  1  -  Pv^v  +  1  <C°’  X)  "v 
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PVVT  1 

V  -  1,  V1  +  1-*V,  V1  nv  -  1  nv'  +  1 


ntA  1  +  A  ]  +  n  A  +  n  A  ,  . 

v  v  -*  v  -  1  v  —  v  -  2  v+2  v  +  2-*v  v  +  1  v  +  l-*v 


n  S9*lm  ,  +  n  ,  SSti,r, 

V  V-*V  -  1  V  -  1  V  -  I—  V 


_stim  stim 

n  S  ,  .  +  n  S  ,  . 

v  v-*v  +  1  v  +  1  v  +  l-*v 


(1) 


The  pumping  and  relaxation  processes  include,  in  order  of  appearance: 
(1)  electron  impact  excitation  of  the  vibrational  states,  (2)  vibration  to 
translation  (VT)  deactivation  processes,  (3)  near  resonant  vibration  to 
vibration  (VVT)  single  quantum  exchange  collisions,  (4)  spontaneous 

radiative  decay,  and  (5)  stimulated  emission  and  absorption  processes. 

VVT  VT  3 

P  and  P  represent  <  vtr>  rates  (in  cm  /s)  for  the  VV  and  VT 

collision  processes.  In  addition  to  the  master  equation  for  the 

vibrational  population  densities,  assumptions  must  be  introduced  for 

calculation  of  the  plasma  parameters  and  radiation  intensities  as  a 

function  of  time  to  theoretically  characterize  the  CO  EDL.  The  sources 

of  the  electrical  excitation  cross  sections,  VT  and  VVT  cross -relaxation 

rates  and  theories,  and  radiative  and  molecular  parameters  used  in  the 

1»  2 

present  model  have  been  discussed  elsewhere. 

The  assumptions  that  are  made  for  determination  of  the  radiation 
intensities  for  an  oscillator  are  the  following.  The  rotational  levels 
are  assumed  to  be  in  thermal  equilibrium  with  the  molecular  kinetic 
temperature,  and  rotational  cross-relaxation  is  assumed  to  be 


sufficiently  fast  that  only  one  rotational  line  can  oscillate  in  any  given 
vibrational  band.  The  gains  for  the  v-*v  -  1  P(J)  transition  are  given 
by 


a(v,J)  =  n  o  (v,  J)  -  n  <7  (v,J), 

v  2 1  v- 1  12 


(2; 


where  cr  (v,  J)  and  <7  (v,  J)  are  the  stimulated  emission  and  absorption 
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cross  sections  defined  by 


*2l(v,J>  = 


J 

S^v)%A10V 

STThT  ^ 


\R  ,/R1n|2  expi-B  J(J-l)/kx| 
1  v, v- 1  1 01  (  v  ) 


<t12(v,J)  = 


S,,V  A  D  J 

v  v  10  v - 1 


BttRT  V 


-JR  ,/Rl2  cxpl-B  J(J!l)/k'Jj  (3) 

1  v,v-l  101  (  v-1  > 
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Losses  such  as  mirror  absorption,  output  coupling,  window  losses,  or 
intracavity  line  selection  losses  are  all  assumed  to  be  uniformly 
distributed  throughout  the  laser  gain  medium.  These  cavity  losses 
are  assumed  to  be  small  enough  that  the  approximation  of  a  spatially 
uniform  intensity  distribution  over  the  length  of  the  active  medium  is 
adequate  for  describing  the  molecular  and  radiative  characteristics 
of  the  medium.  The  optical  resonator  is  assumed  to  be  formed  by  two 
plane  parallel  mirrors  (one  totally  reflecting,  the  other  with  reflectivity 
R)  separated  by  a  distance  L.  From  the  threshold  condition, 

Rexp  [2(  a  -  (3  ) L ]  =  1,  the  gain  coefficient  a.  for  an  oscillating  line 
(for  R~l)  is  given  by  a  -  P  +  (1  -  R)/2L,  where  f3  is  the  loss 


coefficient  for  the  cavity  exclusive  of  the  output  coupling.  The 
percentage  of  the  total  radiative  conversion  of  electrical  energy  that 
can  be  extracted  as  output  (assuming  0  =  constant  for  all  lines)  is 
given  by  the  coupling  efficiency  (T/2)/(  0  L+  T/2),  where  T  is 
the  mirror  transmission.  Radiation  intensities  are  calculated  by 
requiring  that  the  gain  coefficients  o;(v,  J)  of  all  oscillating  transitions 
be  equal  to  the  loss  coefficient  y  at  all  times: 


°2i (v>  J)  nv(t)  “  °i2  (v>  J)  Vi(t)  =  ^  +  T/2L  =  y 


(4) 


or 


Thus,  if  at  time  t  a  set  of  transitions  (v,  J  )  have  just  attained  (< 
continue  to  have)  gain  at  threshold,  their  intensities  lv(t)  must  satisfy 
a  simultaneous  set  of  equations  of  the  form 


<721(v'  V  V*>  -  ai2(v'  Jv>  %.,<*>  = 

-T  [nv(t)  da21(v,  Jv)/dT  -nv  x(t)  dff12(v,  Jy)/dT]  (5) 

The  intensities  I^t)  enter  Eq.  (5)  upon  substitution  of  the  master 
Eq.  (1)  for  nv(t)  , 


n  (t)  =  R  (t)  +  7 

V  V 


hv(v  +  1,  J  ) 

V  +  l 


I  (t) 

V 


hv(v,  J  ) 


(6) 


which  is  here  shown  separated  explicitly  into  the  nonradiative  kinetic 

terms  R-v(t)  and  stimulated  emission  and  absorption  terms.  The 

stimulated  emission  and  absorption  rates  sstim  and  sStim 

V-»  V  - 1  v-l-*v 

that  occur  in  Eq.  (1)  can  be  expressed  in  terms  of  the  cross  sections 
defined  by  Eq.  (3)  and  the  radiation  intensities  IJ  (t)  by  the  relations 


Combination  of  Eq.  (8)  -  (10)  gives  a  relation  for  determining  the 

effective  electron  temperature  T  (t)  as  a  function  of  time. 

e 

If  the  master  Eq.  (1)  is  multiplied  by  the  vibrational  energy  E  and 

v 

summed  over  the  levels  v,  an  equation  for  the  conservation  of  power 
per  unit  volume  is  obtained: 

dEvib<t)/dt  =  we(t)  +  Wvv(t|  +  WV1<t»  +  Wspon<*> 

+  IW.timM  +  Wrot<t>1  (ID 

where  Ey^(t)  is  the  total  stored  vibrational  energy  density,  W  (t)  is 

the  electrical  power  into  vibrational  excitation  defined  by  Eq.  (10), 

W  (t)  ard  W  (t)  are  the  net  powers  converted  to  kinetic  heating  by 

VV  and  VT  collision  processes,  Wgpon(t)  is  fluorescence  power, 

Wstim  is  Power  converted  to  stimulated  radiation  (either  extracted 

or  absorbed),  and  W  (t)  is  power  converted  to  kinetic  heating  by 

rotational  cross -relaxation.  The  last  two  terms  are  shown  grouped 

together,  since  the  origin  of  the  rotational  heating  W  (t)  is  the 

°  rot 

change  of  J  level  that  occurs  with  lasing:  since  stimulated  emission 
and  absorption  of  radiation  occur  on  P(J)  transitions,  the  assumptions 
of  fast  rotational  cross  relaxation  and  rotational  thermal  equilibrium 
with  the  kinetic  temperature  imply  that  the  rotational  energy  change  that 
results  from  (J-l  — *  J)  is  instantaneously  converted  to  heat. 

Explicit  expressions  for  the  terms  occuring  in  Eq.  (11)  are  given  by 

/ 

WstlmW  =  r  2  Mt)  -  y  I(t)  (12) 

V 


WVT(t) 


V,  V'  V  v'  V  v'-l  V,i'-1  —  V  - 1 ,  v' 


V  VT  , 

-n  ,  ii  ,P 

v-1  v'  v-1,  v' —  v,v'-l 


=  £  (E  -  E  n  P  ,  CO,  X 

v  V  V-1  V  V—  v-l' 

X,  v 


-  n  P  (CO,  X)]  N, 

v-1  v-1— V  J  ] 


W  ft) 
rot 


£  (E  -  E  )  [n  SStim  -  n  SSt\m  ] 

V  V  V-1  V  V— v-1  v-1  v-1—  vJ 


w  (t) 

stim 


The  gas  kinetic  temperature  increases  as  a  function  of  time,  due  to  the 
conversion  of  input  power  to  kinetic  heating: 


fNr.n  +  <3/2>  N,J  kT  -  Ww<t>  +  W  ft) 


+  It)  m  WJt) 


(16) 


3-°  approximations  and  scaling  generalizations 


In  order  to  derive  a  useful  set  of  generalized  results  that  can  be 
easily  scaled,  it  is  first  necessary  to  simplify  the  sets  of  equations 
presented  in  the  previous  section  by  introducing  approximations  that 
will  reduce  the  number  of  scale  parameters  required  to  a  reasonable 
minimum  value.  Although  the  additional  approximations  to  be 
discussed  here  are  neither  required  nor  in.-oked  in  the  complete 
computer  analysis,  they  will  be  useful  for  summarizing  the  results 
of  computer  calculations  in  a  coherent  way.  Clearly,  it  will  always  be 
possible  to  characterize  any  particular  class  of  problems  by  a  suitably 
large  set  of  parameters,  so  the  objective  here  will  be  to  show  that  a 
rather  general  class  of  pulsed  CO  oscillators  (containing  CO  and 
monatomic  gases  only)  can  be  described  quite  accurately  in  terms  of 
only  two  basic  scale  parameters  that  define  the  electrical  pumping  and 
cavity  threshold. 


For  a  high  pressure  oscillator  (i.e.,  >20  torr)  operating  at  temperatures 
typically  in  the  range  ~60-300°K,  the  spontaneous  radiation  and  VT 
decay  terms  in  the  master  equation  (1)  can,  to  a  very  good  approximation, 
be  neglected.  Consistent  with  this  approximation,  the  contribution 
kinetic  heating  that  occurs  in  the  temperature  Eq.  (16) 
can  also  be  neglected.  It  will  be  convenient  to  refei  to  the  initial 
values  of  the  temperature  (T  )  and  pressures  (p  ,  n  \ 
in  the  results  of  calculations  to  be  presented  here.  Only  a  small 
percentage  of  electrical  power  is  converted  to  kinetic  heating  (typically 
<  10%  even  after  steady  state  ha?  t  een  attained,  so  (for  reasonably  dilute 
CO  gas  mixtures)  it  is  a  good  approximation  to  neglect  the  temperature 
rise  for  pulse  times  that  are  comparable  to  the  time  required  to  reach 
steady  state.  Thus,  although  the  calculations  correspond  to  the 
situation  for  which  the  gas  number  densities  are  fixed,  the  pressures 


(for  such  cases)  will  change  by  only  a  small  fraction  during  that  time 

scale.  In  any  case,  nothing  is  sacrificed  by  the  choice  ci  the  initial 

values  of  these  quantities  to  form  the  scale  parameters,  and  the  heating 

considerations  will  be  made  more  precise  below.  All  of  the  calculations 

to  be  presented  here  correspond  to  cases  for  which  the  electrical 

excitation  was  kept  constant  as  a  function  of  time  during  the  pulse  -- 

i.  e.  t  -  0,  and  t  =  oo  ,  with  adjustment  of  the  electron  temperature  T  . 
r  f  e 

We  shall  now  show  how  the  coupled  set  of  equations  discussed  in  the  pre¬ 
vious  section  can  be  rewritten  in  a  more  useful  form  by  introducing  the 
following  quantities: 


x 

V 

X 

e 

T 

f 


v  CO 


"e,NCO 


pco' 

1  +  <3/2)ptot/pco 


°Pt/ 


P„<V  .(CO,  CO)  >  <7  (CO,  CO) 
CO  rel _ _ 

£  p  <  v  (CO,  X)  >cr  °pt(CO,  X) 

X  X  rel 


-VVT 


VVT 


=  P  /  (kT  ) 

,v'—  v-1,  v'  +  l  v,  v' —  v-1,  v'+l  mol 


a  , 

V  -a.  V 


=  a  , /{kT  ) 
v-»v'  mol 


ff(v,j)  -  a(v,J)pco/^ 

y  =  >/$ 

w  =  w/p^.0 

i  =  1 * /pco 

E  =  E  /p 

vi.b  vfb  CO 


(17) 
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The  physical  significance  of  these  parameters  is  readily  apparent: 

x^  are  the  relative  populations,  xg  is  a  measure  of  the  degree  of 

ionization,  f  is  r  slated  to  the  number  of  thermal  degrees  of  freedom, 

£  is  the  fractional  percentage  of  optical  line  broadening  caused  by 
-VVT  — 

CO,  P  and  <  v  a  ,  >  are  (proportional  to)  the  VV  and  electrical 

V~*V  -1  -1  - 

excitation  rates  in  units  of  s  torr  ,  a  (v,  J)  in  proportional  to  the 
stimulated  emission  cross  section  with  dependence  on  pressures  removed, 
and  y  ,  W,  I  and  E  are  scaled  threshold  loss,  power,  intensity  and 

•  i 

energy.  It  is  assumed  that  pressures  are  high  er.ough  {  >  20  torr)  that 
optical  lines  are  predominantly  pressure  broadened,  so  that  the  definition 
of  the  parameter  £  is  meaningful. 

Since  the  VT  collision  and  spontaneous  radiation  terms  are  now  to  be 
neglected,  the  master  Eq.  (1)  becomes 


x'(r>=*  £  [*  .  <v* 


v1  v 


>  - 


X  <  V 
V 


o  ,  >] 

v-*v 


+  vfv'  tPv+l,  v’-l-v.v1  xv+l  Vl  +  ■  •  •  I 

-  xv  ?21  (v"J)Tv/hl'v+  Xv-1  Vv’j)  •••  <18> 

where  the  prime  (  )  denotes  differentiation  with  respect  to  T,  and  where 
the  electrical  excitation  term  can  be  written  in  terms  of  the  scaled 
electrical  pumping  power  by  eliminating  the  relative  ionization 

parameter  x  : 


x  =  W  [£  E  (x  <va 
e  e  *-  v,  v  v  v  v 


<  v  a  . 
v 


v>)  r1  (i9) 


Thus,  the  master  equation  expresses  the  evolution  of  the  relative 

populations  x  as  a  function  of  scale  time  T  =  p _ t.  This  is  to  be 

v  CO 
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expected,  since  the  electrical  excitation  is  k'pt  constant,  and  the 

dominant  kinetic  processes  are  (two  body)  V  V  collisions.  Eq.  (18) 

shows  that  the  physics  of  the  molecular  vibration  states  depends 

explicitly  only  upon  the  kinetic  temperature  and  the  scaled  electrical 

excitation  W  =  W  / pl„.  Implicitly,  however,  it  depends  also 
e  e  CO  _ 

the  scaled  threshold  loss  coefficient  Y  for  the  cavity  due  to  the 
oscillation  condition  (4),  and  upon  the  parameter  f  due  to  Lhe  kinetic 
heating  Eq.  (16).  In  terms  of  the  parameters  of  (17),  Eq.  (A)  and 
(16)  become,  respectively 


X  - 


VV,J)  \-l,/,kTmol)  = 


(20) 


T'  =  WH(T)  Tmol/£ 


(21) 


Eq.  (15),  which  defines  the  intensities  of  the  oscillating  transitions, 
becomes 


VV’J'yv  '  VV’J'  X’v-1 


-(WHT  j/OLV  o21(v,J)/dT  -  xv_jd  <712(v,J)/dT]  (22) 


The  power  balance  Eq.  (11)  becomes 


dEib(T)/dT°  5e  +  *VVITI  +  "stim(T»  +  5,rot'T)  ,23) 

It  is  apparent  from  Eq.  (18)  -  (22)  that  the  set  of  parameters 

(T  ,  W  ,  Y  ,  f)  completely  characterizes  the  problem;  all  of  the 
mol  e 

other  (extensive)  input  consists  of  fundamental  physical  constants.  In 
the  following  section,  a  variety  of  numerical  results  based  on  the  complete 
computer  analysis  for  this  model  will  be  presented  and  discussed  to  show 
the  relative  sensitivity  of  these  parameters. 
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4.0 


DISCUSSION  AND  NUMERICAL  RESULTS 


The  computer  calculations  that  have  been  carried  out  correspond  to  a 

fixed  value  f  =  16  (corresponding  to  a  dilute,  10%  CO  mixture).  The 

kinetic  temperatures  T  ^  include  (60,  100,  150,  200,  300)  °K,  the 

scaled  electrical  pumping  rates  (0.  2,  0.  5,  1.0,  2.  0,  5.  0,  10,  20,  50, 

3  2 

100,  200)  W/cm  /torr  ,  and  the  scaled  threshold  loss  coefficients 

(0.1,  0.2,  0.5,  1.0,  3.0,  7.0,  10,  20)  %/cm.  From  the  results  of  any 

given  case  specified  completely  by  the  set  cf  four  parameters 

(T  , ,  W  ,  y  ,  f),  the  general  class  of  problems  that  can  be  scaled 
mol  e 

directly  from  those  results  is  quite  restrictive;  it  allows  only  the 

variation  of  the  CO  partial  pressure,  provided  that  the  electrical 

2 

excitation  W  is  scaled  by  and  that  the  threshold  level  and  gas 

mixture  ratio  are  held  constant.  However,  except  for  the  scaled 

excitation  power  W^,  most  of  these  parameters  do  not  have  a  sensitive 

effect  on  the  power  response  of  the  system  as  a  function  of  time 

(although  detailed  spectral  oredictions  may  have  a  more  critical 

dependence).  Thus,  the  class  of  problems  that  can  be  analyzed  from 

a  given  case  is  often  more  extensive  than  that  just  described.  The 

reasons  for  this  are  the  following:  (1)  the  most  important  parameter, 

the  scaled  excitation  W  ,  is  a  measure  of  the  electrical  pumping  per 

e 

molecule  relative  to  the  rate  of  VV  cross  -  relaxation,  and  thus  specifies 
the  ratio  of  the  two  dominant  kinetic  processes  in  the  CO  EDL.  (2) 

A  pulsed  high  pressure  CO  EDL  typically  operates  in  a  regime  where 
oscillation  is  high  above  threshold,  and  thus,  results  for  optical  power 
extraction  will  usually  be  insensitive  to  the  scaled  cavity  loss 
coefficient  Y  .  Unless  the  ratio  W  /  Y  becomes  very  amall  (to  be 
discussed  more  quantitatively  below),  the  medium  intensity  will  be 
high  above  the  saturation  intensity,  and  the  optical  output  will  be  in¬ 
dependent  of  y  .  (3)  Only  small  percentages  of  electrical  input  power 

are  converted  to  kinetic  heating  (typically  <  10%)  even  after  steady  state 
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has  been  attained,  so  it  is  a  good  approximation  (for  reasonably  dilute 
CO  mixtures)  to  neglect  the  temperature  rise  for  pulse  times  comparable 
to  the  time  required  to  reach  steady  state.  That  is,  the  results  do  not 
depend  critically  upon  f,  providing  it  is  large  (in  which  case,  the 
right  hand  side  of  Eq,  (21)  and  (22)  can  be  neglected).  (4)  Results 
of  calculations  for  systems  operating  with  parameters  in  the  above  ranges 
do  not  show  a  significant  sensitivity  to  dependence  on  T 

mol 

Figure  1-8  show  plots  of  the  quantum  power  efficiency  (for  the  total 
optical  conversion  of  electrical  input)  for  a  variety  of  different  values 

°f  Tmol’  ^  ’  and  We  aS  a  function  of  PCQt.  As  these  families  of 
curves  show,  the  most  important  parameter  for  small  values  of  Y 
is  the  scaled  excitation  power  W  .  Figure  9  displays  a  family  of 
curves  for  which  and  Wg  are  fixed,  with  the  eiiiciency  T)  plotted 

as  a  function  of  PCQt  for  several  values  of  Y  ;  these  show  that  the  optical 
output  is  relatively  insensitive  for  small  values  of  this  parameter.  In 
fact,  Figure  9  suggests  at  least  s emi -quantitatively  that  the  response 
at  100  K  is  approximately  independent  of  Y  if  W  /  Y  >10  W/cm2torr2. 

This  crit^ri-  will  later  be  shown  to  be  related  to  the  scaled  saturation 
intensity  of  the  oscillator. 

It  is  apparent  from  the  set  of  Figures  1-9  that  the  optical  output  efficiency 
displays  a  remarkably  similar  behavior  over  a  wide  range  of  the  parameters 

(We’  Tmol’  *  In  fact'  these  figures  can  be  used  to  show  that,  for  a  given 
temperature  Tq  and  a  fixed  value  of  quantum  efficiency  V /  T)  (where 

T]m  is  the  steady  state  value),  a  loglog  plot  of  p  t  versus  the  scaled 

vU 

power  We  is  approximately  linear.  For  example.  Figures  10-13  show 
such  plots  for  the  turn-on  time  (obtained  by  extrapolating  the  curves 

of  Figures  1-8  to  the  original  and  neglecting  a  small  initial  bump  that  often 
occurs  in  the  calculations)  versus  W  .  The  slope  of  these  curves  is 

? PP roximately  -5/6,  and  is  independent  of  T  ,  and  y  ;  the  variation 

mol 

in  the  magnitude  of  the  curves  is  less  than  a  factor  of  two  within  the 
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range  of  temperatures  chosen.  The  turn-on  time  corresponds  to  the 
point  where  i>=  0;  however,  it  is  apparent  that  the  curves  remain 
similar  as  time  evolves.  Thus,  for  a  fixed  temperature  T  ,  the 

time  t  requii  ed  to  attain  any  specified  value  !»/,_  for  the'optical 
extraction  efficiency  satisfies 


e  = 


,  a  W  5/6 
CO  ?  =  const» 


(24) 


which  implies  that,  within  a  range  of  values  for  7  which  are  suitably 
small,  ^universal"  plot  of  efficiency  versus  the  parameter 

'  =  PCOtWe  can  be  constructed.  Several  examples  of  such  plots 
for  different  temperatures  are  shown  in  Figures  14-17.  |Note  that, 
as  Figures  1-9  show,  the  value  of  the  steady  state  efficiency  r, 
depends  only  slightly  upon  W0,  T^,  and  y  . )  The  significant  of 
the  parameter  e  is  obvious:  if  the  exponent  of  W  in  Eq.  (24) 
were  unity  instead  of  5/6,  «  would  represent  the 'total  electrical 
energy  E  =  W^/p^  (J/Cm3/torr)  which  has  been  deposited  up  to 
time  t.  Therefore,  the  response  of  the  system  is  approximately  a 
function  only  of  the  specific  input  energy  E  (which  is,  of  course,  a 
function  of  time  and  excitation  rate).  The  power  efficiency  >7/  r, 
can  be  integrated  to  give  a  universal  curve  for  the  energy  efficiency 

i  !n„  as  a  function  of  the  parameter  ?,  and  these  results  are 
presented  in  Figures  18-21. 


Among  the  important  conclusions  to  be  drawn  from  the  above  remarks 
are  the  following.  Characteristic  times,  such  as  turn-on  time  or  time 
required  to  reach  steady  state,  satisfy  a  relation  given  by  Eq.  (24). 
Figures  10-17  show  that,  for  the  range  of  temperatures  and  for  the 
reasonably  small  values  of  7  discussed  here,  turn-on  occurs  after  a 
specific  eiectrieai  energy  0.6  .  U  0  .T/liter/torr  has  been  deposited, 
and  steady  state  is  attained  after  an  energy  E  ~  U  ,  .  2.  6  d/liter/torr 

has  been  deposited.  These  numbers  correspond  to  a  specific  pumping 
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~  3  2  — 

rate  W  =1.0  W/cm  /torr  ;  for  other  values  of  W  ,  Eq.  (24)  shows 

6  —1/6  6 
that  these  deposition  energies  increase  as  .  (Note  that 

1.0  J/liter/torr  =  7.  6  is  a  dimensionless  number). 

The  kinetic  temperature  rise  for  pulse  times  comparable  to  the  steady 

state  time  can  be  kept  very  small,  provided  that  the  gas  mixture  is 

reasonably  dilute  (correspondingtto  a  large  value  for  f).  However,  the 

parameter  f  does  not  affect  the  results  critically,  and  is  mainly 

important  for  detailed  predictions  of  spectral  distribution.  To  give 

some  measure  of  the  sensitivity  of  this  parameter,  Eq.  (21)  can  be 

integrated  to  the  steady  state  time.  For  all  of  the  cases  studied,  the 

total  energy  converted  to  kinetic  heating  is  very  small,  typically  <  5% 

of  the  electrical  input.  Thus,  at  steady  state,  the  temperature  rise 

4T  will  be 
ss 


f  4  T  / T  <  .05  E  ~  0.  8 
ss  mol  3  s 


(25) 


Since  the  present  calculations  were  carried  out  for  f  =  16,  this  condition 

predicts  a  total  temperature  rise  limited  to  £  5%.  Clearly,  the  gross 

characteristics  of  the  power  output  are  not  sensitive  to  such  small 

temperature  variations,  so  f  can  be  varied  by  perhaps  as  much  as  a 

factor  of  4  (at  least  for  lower  temperatures  T  ,).  No  detailed 

mol 

investigation  has  been  made  of  the  effects  on  output  spectral  distribution. 


It  is  interesting  to  characterize  the  CO  medium  by  a  saturation  intensity. 
Although  the  full  kinetic  model  contains  oscillation  on  several  lines 
simultaneously,  it  is  convenient  to  neglect  this  feature  for  the  moment, 
and  observe  how  the  gain  saturates  as  a  function  of  the  total  intensity 
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in  the  medium.  From  the  analysis  of  a  simple  two-level  model,  it 
2 

can  be  shown  thac  the  saturated  gain  a  (I)  satisfies 


a  (I)  =  a  (1  +  I/I.)  , 

o  3 


where  the  saturation  intensity  Ig  is  independent  of  the  electrical 

pumping.  We  shall  show  that,  at  steady  state,  the  CO  EDL  follows 

a  simple  dependence  of  this  form,  and  an  estimate  of  the  scaled 

saturation  intensity  as  a  function  of  T  ,  will  be  derived.  In  terms 

mol 

of  scaled  quantities,  Eq.  (26)  becomes,  upon  setting  O' (I)  =  y , 


V I  y  I/I  =  a  . 

s  o 


At  steady  state,  y  I  =  rj  W„,  and  thus, 

00  “ 


y  +  fl  w  /I  =  a 

00  e  s  o 


This  equation  shows  that,  for  a  fixed  value  of  W  ,  the  steady  state 
power  efficiency  and  the  cavity  threshold  V  are  linearly 

related.  That  is,  if  ^^(O)  i3  the  limiting  value  of  the  efficiency 
as  y-»0,  then 


7  *  -  ( V  VoiKA 


Figure  22  shows  tj  plotted  versus  7  for  T  =  100°K,  for  two 

'Z.  mol 

different  values  of  W^,  and  as  Eq.  (29)  shows,  the  slope  |  d  y  /drj^  |  of  thes 
curves  is  proportional  to  Wg.  Furthermore,  the  proportionality  constant 
determines  the  scaled  saturation  intensity,  which  is  determined  from 
Figure  22  to  be  Tg~0.8  W/cm2/torr2.  Similar  plots  for  a  fixed  value 
of  We  and  variable  temperature  are  shown  in  Figure  23.  All  of 

these  curves  are  linear  to  a  very  good  approximation,  and  the  saturation 
intensities  obtained  from  their  slopes  are  summarized  in  Table  I. 


During  the  transient  evolution,  the  relation  between  the  quantum 

etficiency  and  the  cavity  threshold  dees  not  seem  to  follow  a  linear 
relationship,  as  the  results  of  Figure  9  indicate.  This  case  is  shown 
in  Figure  24,  in  which  the  efficiency  curves  for  three  values  of  y 
have  been  replotted,  with  efficiencies  at  specific  values  of  p  t  plotted 
versus  V  alongside  for  comparison. 

Table  I.  Scaled  Saturation  Intensity 


Tmo)  <°K> 


I  (W/cm^/torr^) 
s 


60 

100 

150 

200 

300 

0.  6 

0.  8 

1.6 

2.  5 

6.  4 

Small  signal  gain  characteristics  can  be  obtained  from  Eq.  (28)  with 
7  =  0*  For  values  of  pumping  and  cavity  threshold  considered  here, 
it  is  apparent  that  the  small  signal  gain  is  proportional  to  We, 

•<*n  =  (0)  W  /I 

°  e  s 

For  efficient  operation  of  the  osc  llator  (high  above  threshold),  the 
necessary  condition  is  aQ»T  »  a...d  therefore  (since  r]  —  1), 


W  /y  »  I  . 
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Quantum  efficiency  as  a  function  of  scale  time 
several  values  of  scaled  electrical  pumping  power 
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Quantum  efficiency  as  a  function  o  f  scale  time  Pqq*  for 
several  values  of  scaled  electrical  pumping  power. 
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Figure  14.  Quantum  powei  efficiency  as  a  function 
of  the  parameter  £  =  p^tW5'6.  Note 
that  e  is  approximately  energy  deposited 
in  units  of  J/liter /torr. 
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Figure  15.  Quantum  power  efficiency  as  a  function 
of  parameter  e  =  p  fW5'6.  Note 
that  e  is  approximately  Energy  deposited 
in  units  of  J /liter /torr. 
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Figure  16.  Quantum  power  efficiency  as  a  function 
of  parameter  e  =  P^Q^e  -^ote  t^iat 

e  is  approximately  energy  deposited  in 
units  of  J /liter /torr. 


Quantum  Power  Conversion  Efficiency  T)(  t)/7»a> 
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1.0  1.5  2.0  2.5 

Pco1^76  tms* t°rr(W/cm^/torr2)^^] 

Figure  17.  Quantum  power  efficiency  as  a  function 
of  parameter  e  =  '  Note  that 

e  is  approximately  energy  deposited  in 
units  of  J/liter/torr. 
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Quantum  Energy  Conversion  Efficiency 


(ms*torr*  (W/cm  /torr  )  ) 


Figure  19. 


Quantum  energy  efficiency  as  a  function  of  paramete 


Note  that  e  is  approximately  energy 


e  =  p 


depos:Fed  in  units  of  J /liter /torr. 


20.  Quantum  energy  efficiency  as  a  function  of  parameter 
e  =  PCQtW  .  Note  that  e  is  approximately  energy 
Qeposited  in  units  of  J/liter  torr. 
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Quantum  Energy  Conversion  Efficiency  £(  t)/T) 


1.0  2. 0  3.0  4.0  3.0 

p  t  W  (  ms.torr.  (W/cra^/torr^)^^) 

w  6 

Figure  21.  Quantum  energy  efficiency  as  a  function  of  parameter 
e  =  p  tW5/  .  Note  that  e  is  approximately  energy 
deposited  m  units  of  J/liter /torr. 
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Scaled  ThreshholcI  Loss  Coefficient 


Saturation  characteristics  of  a  CO  EDL  after  attainment  of 
steady  state;  the  slope  of  the  curve  is  a  measure  of  W  /I  . 
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Steady  State  Quantum  Power  Conversion  Efficiency 


Scaled  Threshhold  Loss  Coefficient  y  (%/cra) 

Figure  23.  Saturation  characteristics  of  a  CO  EDL  after  attainment  of  steady 
state;  the  slope  of  the  curve  is  a  measure  of  W  /I 
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igure  24.  Saturation  characteristics  of  a  CO  EDL  investigated  as  a  function  of  time.  Curves  on  left 
are  plots  of  quantum  power  efficiency  n(t)  for  three  values  of  V,  and  on  right,  q(t)  is 
plotted  versus  V  for  different  values  of  p  t. 
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